Description of the Oxygen Order evolution and its relation to the Superconducting 

Transition in La2Cu04^y 



E. V. L de Mello 

Instituto de Fisica, Universidade Federal Fluminense, 
Niter 6i, RJ 24210-340, Brazil 

(Dated: February 21, 2012) 

The segregation of oxygen in the high critical temperature cuprate superconductor La2CuOA+y 
has been systematically studied along the years. In a recent set of experiments, Poccia et al related, 
for the first time, time ordering (t) of oxygen interstitials with the corresponding superconduct- 
ing transition temperature Tc(t). We develop a phenomenological description of the time ordering 
forming pattern domains and show how it may affect the superconducting interaction. The super- 
conducting self- consistent calculations in a system with electronic granular structure of varying hole 
doping yields also different local d-wave amplitudes. These amplitudes are of the order of mag- 
nitude of scanning tunneling microscopy measurements and they vanish at T*(t) > Tc(t). Then, 
calculations with Josephson coupling among the isolated charge domains reveal that the supercon- 
ducting interaction is likely to be scaled by the local free energy and capture the details of Tc(t). 
The accurately reproduction of these apparently disconnected phenomena establishes routes to the 
important physical mechanisms involved in the connection between sample production and on the 
origin of the superconductivity of cuprates. 

PACS numbers: 74.20.Mn, 74.25.Dw, 74.62.En, 74.81.-g 



I. INTRODUCTION 

There are considerable evidences that the tendency to- 
ward phase separation is an universal feature of some ox- 
ides like cuprate superconductors and manganites[l|, Q. 
The presence of hole-rich and hole poor phases were de- 
tected in oxygenated La2Cu04-\-y almost immediately 
after the discovery of the superconductivity in these 
compounds [3] and by subsequent works [^-^Sj. Further 
experiments have observed evidences of complexity and 
electronic disorder in many other materials [l|, H, [9[ . 

To describe this phenomenon in the context of cuprates 
some theories producing phase separation have been sug- 
gested, mainly based on doped Mott-Hubbard insula- 
tors. To cite just a few, some relies on Fermi-surface 
nesting which leads to a reduced density of states, or a 
gap, at the Fermi energvfiol fllj. Others uses a competi- 
tion between the tendency of an antiferromagnetic insu- 
lator to expel doped holes and the long range part of the 
Coulomb interaction to explain the formation of charge 
ordered phases [H, [l^ . Another approach suggested that 
the pseudogap line T* in the onset of a 1st order transi- 
tion, with the development of carriers of two types, frus- 
trated by the electro-neutrality condition in the presence 
of rigidly embedded dopant ions [3, [HI . To deal with this 
problem some authors considered the presence of impu- 
rities or inhomogeneities in the local electronic distribu- 
tion as a starting point to derive the superconducting 
properties [l6l, fl7|. Some others approaches used the in- 
fluence of a mesoscopic phase separation in the appear- 
ance of superconductivity to derive the general properties 
of cuprates 

These theories describe the observed features of some 
properties, however they fail to predict some basic ex- 



perimental results, specially those related with some 
form of real space inhomogeneities. Specifically, the re- 
cent combined experiment relating the time evolution 
(t) of the domain growth of oxygen interstitials (i-0) 
in La2Cu04-^y with the subsequent measurement of the 
superconducting transition temperature Tc(t) [20[ brings 
new possibilities that require new approaches. 

In this paper we describe the two parts of the Poccia 
et al experiment [20] and try to provide an explanation to 
the raise of Tc with the oxygen ordering. We rely on a 
description of the phase separation in cuprates based on 
the time-dependent Ginzburg-Landau (TDGL) or Cahn- 
Hilliard (CH) equation introduced earlier [TqI |21|. It is 
well known that these equations display solutions with 
two type of domains where the free energy assumes some 
minimum value, for high and low density. These valleys 
are surrounded by steep boundaries where the charges 
can get trapped, loosing part of their kinetic energy what 
enhances the mechanism of pair formation [22l-[2^ . The 
free energy barrier between the two phases varies with 
the time what connects the time of oxygen ordering with 
the variations on Tc(t). Scaling the superconducting in- 
teraction with the local changes of the free energy is one 
of our most interesting finding. Another new point is 
that in this granular-like system the resistivity transition 
temperature Tc occurs when the Josephson energy Ej 
among the grains is equal to Kb^c[23|, [23 • These new 
ideas are endorsed by the close agreement with the data. 

It has been determined experimentally that the su- 
perconducting state occurs only at some values of y in 
La2C2i04+y [HI,!!! it is enhanced when the excess 

oxygen dopants form a three-dimensional ordered super- 
lattice within the interstitials regions of the crvstal(8[[25j. 
Systems with poor i-0 ordering may have two types of 
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FIG. 1. (color on line) Plots taken from Poccia et al[20]. On 
the left panel the presence of the reflected X-ray diffraction 
peaks associated with the oxygen order which, after a thresh- 
old, increases rapidly and enter in a low variation regime 
after a few hours. After being illuminated the samples are 
quenched to low temperatures and, on the right panel it is 
plotted the corresponding measured values of Tc{t). 



superconducting states with lower Tc's than the one with 
more i-0 ordering[25|. On the other hand, since the su- 
perconductivity occurs mainly in the Cu02 planes, it is 
clear that there is a connection between the i-0 oxigen 
order and the charges in the planes. This connection has 
also been observed in other cuprates where the charges 
tends to aglomerate near the dopants Consequenly, 
it is an open question whether the electronic disorder is 
driven by the lower free energy of undoped antiferromag- 
netic (AF) regions [22| (intrinsic) or by the out of plane 
dopants (extrinsic origin) [26], like the i-0 ordering[8, 2^. 
Notice that the i-0 disorder temperature (330K) is quite 
similar to Neel temperature and perhaps there is a con- 
nection between the two phenomena [22|. In either case, 
being intrinsic or extrinsic, it is possible to assume that 
there is an one to one correspondence between the two 
phenomena. This one to one correspondence is a crucial 
ingredient of our work and will be explored in detail here. 
Consequently the i-0 ordering time evolution is assumed 
to be described together with the planar electronic phase 
separation (EPS) by the CH equation [li, [2l|-[23|. 



will turn our attention the hole dynamics in the Cu — O 
planes. Using the TDGL or CH approach, the appropri- 
ate order parameter is the normalized difference between 
the local p(z, t) and the average electronic doping level p 
u{p^i^t) = {p{i^ t) —p)/p. Here, in order to compare with 
the Poccia's experiment, we assume the optimal doping 
p = 0.16. Clearly u{i^t) = corresponds to the homoge- 
neous system above Tp^, and u{i^t) = ±1 correspond to 
the extreme cases of full phase separation that probably 
is not the case for the superconductor La2Cu04-^y^ be- 
cause the i-0 mobility occurs only in a small temperature 
interval below the phase ordering temperature. In other 
similar system like La2-xSrxCuO/^ the degree of phase 
separation seems to be larger since AF domains (with al- 
most zero local doping) are observed even in overdoped 
compounds [2 7]. Consequently, in the i-0 experiment it 
is likely that the system reaches an intermediated struc- 
ture with two phases with densities around the average 
doping level p[25]. 

Then, in terms of the order parameter u{p^i^t), the 
Ginzburg-Landau (GL) free energy functional is the usual 
power expansion. 
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Where 

B'^u^/A + A^(T) = a{Tps{p) - T), a and 5 are con- 
stants, e gives the size of the boundaries between the 
low and high density phases [HI, |2l| . The CH equation 
can be written [28j in the form of a continuity equation 
of the local density of free energy /, dfU = —V.J, with 
the current J = M\/{Sf/5u), where M is the mobility 
or the charge transport coefficient that sets the phase 
separation time scale. Therefore, 



du 
'dt 



-MV\e^V^u - A\T)u + B^w^). (2) 



II. THE PHASE SEPARATION DESCRIPTION 

Normally the i-0 order starts below the phase separa- 
tion transition temperature Tps ^ 330K. In the Poccia 
et al^ experiment the system is heated to a disordered 
phase above 370K and quenched to 200-300K where the 
i-0 have low mobility and the sample exhibits no i-0 or- 
dering and a poor superconducting state. However, by 
illuminating such a system at this range (200-300K) of 
temperature, they observed nucleation and growth of the 
ordered i-0 domains with the time of irradiation t. Here 
we assume that this process is accompanied also by the 
nucleation of charge domains in the planes. A similar 
process where the holes in the planes follow the dopant 
atoms and segregate in inhomogeneous domains was ob- 
served in YBCO|26j. 

Based on these facts, and that the superconducting 
transition is accomplished by the electronic charges, we 



The parameter e determines the planar size of the 
boundaries and A{T)/B the values of the order pa- 
rameter at the minima near the transition temperature. 
If A(T) is zero (above the phase separation temper- 
ature Tps) there is only one solution for the free en- 
ergy and for the non-vanishing case there are two so- 
lution corresponding to the two phases (high and low 
densities) (iil . In Fig.(j2j), we show same typical simu- 
lations of the density map with the two (hole-rich and 
hole-poor) phases given by different colors. Varying the 
values of the ratio A(T)/B and the symmetry of orig- 
inal seed we can obtain these different patterns. In 
some previous works [l9[ [22|, [29| we have tried to match 
these patterns with some experimental evidences, like the 
patchwork forms of the scanning tunne ling microscopy 
(STM) local density of states (LDOS)[30|43il or the stripe 
like patterns as measured in some systems by neutron 
diffractionfo^. In theory the ratio A{T)/B is zero at Tps 
a should increases as the temperature decreases but since 
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FIG. 2. Different patterns derived from Eq.(|2| with different 
parameters and different seeds. In all cases the local order 
parameter evolves from zero to zLA/B. The size between the 
walls separating the two phases is scaled by e. We choose 
to work with the pattern of the buttom-right panel because 
it ressembles the phase separation picture of Ref . ( [25l] ) amd 
Ref.([3l). 



the mobility also decreases with temperature we take a 
constant ratio A/B in the simulations. For instance, in 
the case of the La2Cu04-\-y experiment the i-0 ordered 
phase develops only in the small interval between 330K 
and 300K[20||. 

The simulations to the i-0 experiment were made with 
£ = 0.08 and A{T)/B = 1 that gives the patterns shown 
in the button-right of Fi g. ^ because it resembles the 
STM patterns on LASC0[3j. As in the case of YBCO 
where the holes seem to follow the dopant atoms [26|, we 
assume that the planar charge follows the i-0 and a sim- 
ilar phase separation occurs within the planes. 

In the Poccia et al experiment pof the system is heated 
to T = 370K above the mobility temperature Tm = 330K 
and quenched to low temperature. Such system does not 
present ordered i-0 peaks and has poor superconducting 
order. However, by illuminating this disordered system 
with X-rays, they observed, after a time threshold to, 
the nucleation and growth of ordered domains which is 
accompanied of the recovery of a robust high Tc state. 

Therefore we simulate the appearance of the domains 
by the formation an EPS structure. In the simulations 
on a 105 x 105 sites, we note that below t=250ts (ts=time 
steps) the solutions of Eq.([2j) have an uniform density but 
at t=265ts a disorder with a regular checkerboard pattern 
sets in. This solution with high and low densities evolves 
up to 2000ts when another domain pattern, more irreg- 
ular, sets in after a transient time as it is shown on the 
top-right panel in Fig.®. This pattern is consequence of 
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FIG. 3. (color online) The time evolution of the density 
view map simulation of the inhomogeneous charge density 
with average p = 0.16 on an array of 105 x 105 sites. The 
vertical scale gives the color code for the local densities (red 
is high and blue is low density). Top panels are for t = 2000 
end t = 4000ts. Lower panels are for the t = 20.000 and 
t = 200.000ts. In top middle the local density p{i) histogram 
for the t = 2000ts = 0.1 /i because the stable pattern start 
to develop after this transient time. In the lower panel the 
histogram for t = 200.000ts = lOh. The estimated FWHM of 
the charge dispersion A of each low and high density phases 
are marked in the plots and are to be compared with the i-0 
ordering Bragg peaks experimental dispersion values. 



more stable solution which will dominate at larger times. 
Consequently, we take this solution to match the order- 
ing promoted by x-Ray illumination. Thus this stable 
solution, starting at the transient t=2000ts should corre- 
spond to the i-0 ordering phase where the Bragg peaks 
sets in. 

At 4000ts, the checkerboard order rests only on less 
than 1/3 of the system and the more stable irregular 
granular order dominates as shown in Fig. (j3j. This might 
be the situations described by Fratini et al[25|, a system 
with two phases and different TcS. Above 6000ts the 
granular pattern dominates and from 20000ts (button- 
left of Fig.®) to 200000ts (button-right of Fig.©) the 
grains grow very slowly. To compare with the exper- 
iments, we take the to = 2000^5 or O.lh, that is the 
onset of the stable phase, that one that will remain at 
larger times, as the onset or threshold time of the raise of 
the diffraction peaks associated with the ordered phase. 
This is because at to = 2000ts that the stable phase 
sets in and, in the experiment, at to is when the i-0 
ordering develops and remains at larger times. This con- 
nects in a phenomenological way the oxygen ordering at 
larger times with the raise and grow of the stable do- 
mains shown in Fig.(|3]). In the CH approach, another 
way to follow a phase separation process is through the 
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FIG. 4. (color online) The FWHM of the i-0 X-rays 
diffracted dispersion peaks [S^] and the values of the local 
high and low local densities histograms dispersion A(t) cal- 
culated from 6 different times density maps dispersion (blue 
stars) as shown in Fig.JSl 



local densities histograms. In Fig. ([3]) we show this pos- 
sibility for two selected times and as expected, the low 
and high local density dispersion A{t) decreases as time 
increases. This is equivalent to the phase separation into 
two distinct phases, evolving from just one average dop- 
ing phase. We compare directly the two dispersions, the 
calculated and the experimental data in Fig.(j4]), without 
any adjustable parameter. 



for the case of t=10h. It is also shown in the left panel 
the values of VcLiht) along 25 sites in a straight line in 
the middle of the view map of Fig.(j5j), for four different 
times of phase separation. In this way, we can visualize 
the regions where the charges get trapped and the inset 
shows the time variations of the barrier walls that, by 
assumption, scale the superconducting inter action [2^. 
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FIG. 5. (color online) A 3D view map of the potential 
VGL{i,t = lOh) on the right panel. On the left, the values in 
row on 25 sites i to visualize the free energy potential wells 
with walls defined as Vgb{t) in the text. The values of Vgb{t) 
are shown here at t = 0.3, 3/i, 5h and lOh and are very similar. 
In the inset we blow up the peaks and we see clearly an steady 
increase between t = 0.3 and t = 5h and the saturation above 
t = 5h. Notice that the button of the energy curves, corre- 
sponding to the grains, remain fixed. The relative variations 
on Vgb{t) are very small but have an important role in the 
BdG calculations of Ad{t) and on the Tc{t). 



III. THE LOCAL FREE ENERGY AND THE 
PAIR FORMATION 

The phase separation process described by the CH dif- 
ferential equation occurs due to the minimization of the 
free energy given by Eq.([T]). Solving the CH equation, the 
potential VcLiht) can be calculated directly from free 
energy (Eq.([T])) through the order parameter u{i^t). We 
performed a detailed study of this potential and found 
that it takes it minimum value at the domains and that 
is why the domains are formed. We notice that, as the 
time evolves, the minimum value remains practically un- 
changed but the height of the barriers increases slowly. 

On of the main results of Poccia et al^^ is that the 
samples with more time of exposition to X-rays devel- 
ops i-0 ordering regions and have larger Tc. This is 
a strong indication that the formation of i-0 ordering 
domains is associated also with the formation of charge 
domains in the planes where the superconducting proper- 
ties develops. Consequently, the trapping of the charges 
in inhomogeneous domains may be connected with the 
superconducting properties. To explore this possibility 
we study in detail the energy maps of the VcLiht) for 
different times t 

To start this analysis we plot a 3D view map in Fig.© 



The calculations shown in the inset of Fig.© indicates 
the possibility to connect the phase separation time with 
the height of the VcLiht) between the domains walls, 
which we define as Vgb{t). Such free energy barrier is re- 
lated with probability of the planar charges be trapped 
in the domains as shown in a more clear way in Fig. ([6]). 
The effect of Vgb on the charges is to attract them to the 
grains, loosing part of their kinetic energies. This loss of 
kinetic energy in the presence of a two-body attractive 
potential favor the possibility of Cooper pair formation 
and the domain walls acts as a catalyst to the supercon- 
ducting state [23I, [IJ. At this stage of our phenomeno- 
logical approach we are not able to decide if Vgb is the 
origin of the pair potential or if just strength some other 
fundamental interaction. The small changes on the size 
of Vgb{t) with the time of phase separation t, as shown if 
Fig.©, may be the reason of the changes in the super- 
conducting amplitudes that reproduces the experimental 
values Tc{t). Consequently, we assume that Vgbit) scales 
the pair interaction used in the superconducting calcula- 
tions of next section. 

The definition and derivation of the values Vgb{t) is 
essential to our calculations since its variation will be 
used in the Bogoliubov-deGennes calculations below to 
connect the time t of x-ray illumination with the super- 
conducting critical temperature Tc{t). Therefore, we plot 
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in the Fig.(j6j) how it is defined from a typical potential 
energy simulation at a given time. 




FIG. 6. (color online) The definition of the values of Vgb, 
related with the height of the potential barrier among the low 
and high energy grains is shown in the top inset. This plot 
is made following Vohii^t — 6200) along the white straight 
line in the middle of the figure. As the time (of irradiation) 
increases, so does the height of the barriers (see Fig.©) and 
the probability of the charges become confined also increases. 



where fn is the Fermi function. We stop the self- 
consistent calculations only when all p(x^) = p{i) con- 
verges to the appropriate time CH density map like those 
shown in Fig.(|3j). 

Typical solutions can be visualized in Fig. ([7]) where 
we show the local density p{i) on a square of 28x28 sites 
where the BdG calculations were made and the 3D map 
of the low temperature /S^dii^T ^ Oi^). The pairing po- 
tential Vgi){t) ^ tij = O.lSeV is parametrized to match 
the average local density of states (LDOS) gaps measured 
by low temperature STM on the optimal doping p = 0.16 
LASCO [Ili, namely, Ad{T ^ 0) ^ 7 - 12meV. Starting 
with Vgb{t = 2000ts), analyzing the results of Vgb{t) like 
those shown in Fig.(j5j), we were able to connect the time 
of x-ray exposure to the superconducting amplitudes. 
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IV. THE SUPERCONDUCTING 
CALCULATIONS 

The CH simulations yield density maps p{i^t) on a 
square lattice like that of Fig.(|3]). They are used as the 
initial input and it is maintained fixed throughout the 
self-consistent Bogoliubov-deGennes (BdG) calculations. 
We use nearest neighbor hopping tij = 0.15eV and next 
nearest neighbor hopping of ti/tij = —0.27 taken from 
hole doped experimental dispersion relatio ns [331 . For 
completeness, the BdG equations areplol [22I. 



(3) 



These equations, defined in detail in Refs.pJl, |3^J, are 
solved self-consistently. Un, Vn and En > are respec- 
tively the eigenvectors and eigenvalues. The d-wave pair- 
ing amplitudes are given by 





Ad(xO 



[Un{^i)v*{Xi + (^) + <(Xi)iin(Xi 



+(5)] tanh 



En 



2kBT' 

and the local inhomogeneous hole density is given by 



(4) 



p(x,) = 1 - 2^[|u„(x,)|2/„ + |«„(x.)|2(l - /„)], (5) 



FIG. 7. (color online) Left panel shows a typical density view 
map used here on a 105x105 sites and the 28x28 square where 
the local superconducting amplitudes Ad{i,t) are calculated 
by the BdG approach at t = GOOOts = 0.3h, as an example. 

With these values of Ad{i^T ^ 0), we find that 
A(i{i^T) ^ at a single temperature T* ^ 96K which 
is much larger than typical values of Tc{p = 0.16) 
for La2Cu04-^y^ but in agreement with the pseudogap 
phase measured by STM[3i|, or the Nernst signal on 
LSCO[29|, [35|. The variation of the superconducting am- 
plitudes with temperature is shown in Fig.®. Since this 
value of 96K is much larger than the values of for 
optimal La2Cu04-\-y^ we notice that the EPS transition, 
with the free energy walls and wells, makes the structure 
of the system similar to granular superconductors (36j . 

In this way, it is likely that the superconducting transi- 
tion occurs in two steps [isl,!!^: first by intra-grain super- 
conductivity and than by Josephson couphng with phase 
locking at a lower temperature. This approach provides a 
clear interpretation to the pseudogap phase of cuprates. 

To obtain the experimental values of Tc(p), we use 
an approach similar to a system composed of gran- 
ular superconductors (36) . Our main proposal is that 
the superconducting transition occurs in two steps as T 
decreases [23, 24]: first by intra-grain superconductivity 
and than by Josephson coupling with phase locking at 
low temperatures. 

These two completely different calculations, yielding 
two different energy scales, are motivated by the two en- 
ergy scales found in most cuprates mentioned above [? ? 
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FIG. 8. The values of Ad(i, T) x T at 4 different points "i", 
for the case t = 6000ts = 0.3h. The values at low temperature 
agree with the values measured by the LDOS from the STM 
data on LSCO. Notice that they vanish near T* = 96K much 
above the measured Tc ~ 35K. 



] , and also in the two different regimes of the fluctuation 
magnetoconductivity in a YBa2Cus07 single crystal (37j. 
These systematic measurements detected an effectively 
two-dimensional (2D) regime far above that can be 
interpreted in our approach by the isolated superconduct- 
ing grains in the CuO without phase coherence. Decreas- 
ing the temperature towards Tc a crossover to a three- 
dimensional (3D) Gaussian regime sets in as the different 
superconducting regions or grain become connected and 
increase coherence. Very close to Tc a critical regime 
characteristic of a 3D XY university class is measured 
which is consistent with a granular superconductors in 
which each superconducting grain develops its own phase 
^ that may oscillate but become locked at Tc, yielding 
long range order. 

Therefore, due to the formation of the local or intra- 
grain superconducting amplitudes at low temperatures, 
the system is regarded as an array of Josephson junc- 
tion as shown schematically in Fig.(|9]). The proper- 
ties of Josephson tunnel junctions between d-wave su- 
perconductors have been studied by Bruder et aljssj. 
They calculate the tunnel matrix elements in second- 
order perturbation theory for two superconductors (1 
and 2) with superconducting amplitude 1/2(^5 T, <l>) = 
A{i^T)cos[2{^ — ^1/2)]- We use their results to an array 
of junction as schematically shown in Fig. ([9]), where the 
local amplitudes Ad{i,T) are represented by the size of 
orbitals and the phase angles (j) have been drawn around 
x-direction. <l> is the polar angle of the flrst Brillouin 
zone. 

To calculate the Josephson coupling energy that is con- 
nected with the onset of phase coherence, we recall the 
work of Bruder et al[38| on d-wave superconductors junc- 
tions. They found that the tunneling current behaves in 
a similar fashion of s-wave superconductors junction and 
the leading behavior is determined by tunneling from a 



FIG. 9. Schematic figure representing the array of Josephson 
junctions on a Cu02 plane. The local dependent amplitudes 
A(p, i, T) are represented by the orbitals and the phase angles 
(j) oscillates around the crystal a-axis direction due to the 
symmetry of the d-wave on the Cu02 plane. 



gap node in one side of a junction into the effective gap in 
the other side. Consequently, as a flrst approximation to 
the Josephson coupling energy Ej^ we adapt the theory 
of s-wave granular superconductors [39] to the A^^(j9, T) 
in the grains. 

Where A»^(T,t) = Ad(T,i,t)/Ar, since Ad(T,i,t) 
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FIG. 10. (color online) Left panel, the thermal energy ksT 
and the Josephson coupling among superconducting grains 
Ej{t,T) (as given by Eq.®) for some given time of irradia- 
tion t as function of T. The intersections give the calculated 
results for Tc(t) plotted with inverted triangles on the top of 
the experimental results in the right panel. 



varies with the position i in the inhomogeneous sys- 
tems, as shown in the right panel of Fig. ([7]). Rn is the 
normal resistance of the La2Cu04-^y compound, which 
we assume to be independent of the time t as inferred 
from the data of Poccia et al[20]. It is also propor- 
tional to the planar resistivity pab measurements [40] on 
the La2-pSrpCu02 series. In Fig. ([10]), the Josephson 
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coupling Ej{p^T) is plotted together with the thermal 
energy KbT. The intersections yield Tcit). 

As discussed in connection with the free energy of 
Eq.([T]) and Fig. ([5]), the relative values of the grain bound- 
ary potential wall Vghit) are easy to estimate, but we do 
not know their absolute values. Consequently we use 
Vnh (t = 0.2/1) that matches the 1st point of Poccia et 
aljiol, namely, Tc{t = 0.2h) = 33.6K. Ah the others 
Tc{t) results for larger t, shown in Fig. ([10]), follow by 
the small but steady variations in Vgb{t) shown in the in- 
set of Fig.(j5j). In other words, we use Vgb{t = 0.2h) as a 
parameter to calculate Tc{t = 0.2h) and all the others 7 
points shown in the right panel of Fig. ([10]) are parameter 
free, that is, they are calculated taking the relative vari- 
ations in the free energy barriers as shown in the inset 
of Fig. ([5]). The steep increase of Tc{t) during the consol- 
idation of the stable phase between t ^ 0.3 — 5h agrees 
well with the data as shown in Fig. ([TQ|) . indicating that 
our triple calculations (CH simulation of BdG cal- 
culation of local Arf(z, t) and the Josephson estimation of 
Tc{t)) provides a new way to describe the assembled of 
superconducting properties in La2Cu04-\-y. 

V. CONCLUSIONS 

In conclusion, we have provided a complete new de- 
scription of a set of complex experiments relating the 
time of irradiation that booster i-0 in ordering domains. 



with the superconducting temperature Tc{t). The varia- 
tion of Tc{t) shows that this out of plane phase separation 
induces also changes or phase separation of the charges 
into the Cu — O planes. Consequently, we use time de- 
pendent CH approach to describe the associated charge 
segregation in the planes. This method allows also to 
follow the local variations in the free energy during the 
phase separation process. Analysis of the local variation 
of the GL potential in the charge inhomogeneous pro- 
file shows that the charge gets trapped into the high and 
low density domains. The confinement of the holes in 
these domains or puddles reinforce the possibility of lo- 
cal Cooper pair formation in these isolated domains as 
in a granular superconductor where different grains may 
be connected through Josephson coupling at low tem- 
peratures. Assuming that the superconducting interac- 
tion in the puddles may be strengthen or scaled by the 
height of the free energy barriers Vghit)., we were able to 
connect the superconducting Tcit) through the BdG ap- 
proach on a charge disordered profile in the planes with 
the time of X-ray illumination. This procedure allows us 
to connect for the first time the degree of phase separa- 
tion or the degree of inhomogeneities generated at high 
temperatures with the low temperature superconducting 
properties. These findings reveal a virtually unexplored 
line of research on how systematic variations of sample 
preparation may affect the superconductors properties. 
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